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HIGH SUBSONIC FLOW TFSTS OF A PARALLEL PIPE FOLLOWED 
BY A LARGE AREA RATIO DIFFUSER 


By 

P.S. Barna^ 

SUMMARY 

Experiments were performed on a pilot model duct system in order to explore 
its aerodynamic characteristics. The model was scaled from a design projected 
for the high speed operation mode of the Aircraft Noise Reduction Laboratory 
located at NASA Langley Research Center. The test results show that the model 
performed satisfactorily and therefore the projected design will most likely 
meet the specifications. The experiments were conducted in the Aerodynamic 
Laboratory of Old Dominion University. 
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LIST OF S'/MBOLS 


A = sectional area, 

a = throat area of Venturi, 

0 

Cp = diffuser pressure recovery factor 

c, = coefficient of discharge 
d 

D = diameter of the parallel pipe, m 
7 = average friction factor along parallel pipe 

g = gravitational acceleration, m/sec^ 

K = specific heat ratio 

AL = length of pipe section under consideration 
M = Mach number measured along parallel pipe and diffuser 
p = absolute pressure measured along parallel pipe and diffuser, kg/m 
Ap = pressure difference in Venturi, kgf/m^ 

Q = volumetric flow rate, m^/sec 

R = gas constant, joule/kg ®R 

T = stagnation temperature, “R 

V = velocity, m/sec 

W = flow rate, kg/sec 

X = distance measured along parallel pipe and diffuser 
y = distance across the flow, m 

Y = specific weight, kgf/m^ 

9 = static temperature 

p = fluid density, kg/m^ 

p = fluid density at isentropic stagnation conditions, kg/m^ 
o 

Subscripts 

1 = inlet to a section 


2 


outlet 


INTRODUCTION 


The experiments presented in this report were planned to explore the flow 
behavior of an air duct system scaled down to approximately a one-tenth* size 
pilot model of the prototype. The prototype itself is a duct system projected for 
the Acoustic Noise Reduction Laboratory at NASA Langley Research Center, The three 
main components of this system are an air intake, a parallel pipe, and a diffuser. 
The flow entering the •’..‘r intake is to be accelerated in transit to a Mach number 
of about 0.5 and is to be further accelerated along the parallel pipe to unity 
(sonic condition) at pipe exit where the flow is permitted to enter the diffuser 
for pressure recovery. A large centrifugal blower is projected to perform the 
necessary flow induction, and it is anticipated that the model operation will 
furnish satisfactory data for estimating the performance of the prototype. 

Preliminary studies indicate that there arises no particular problem with the 
intake and a sensible design will satisfy the need for a uniform distribution of 
flow tc the pai'allel pipe. The requirement of attaining a Mach number 0.7 along 
the pipe, and preferably being able to increase it to unity at pipe exit, calls 
for attention considering the variation of friction factor along a parallel pipe 
anticipated under hij.My accelerating flow conditions. Furthermore, the uncon- 
ventional design of the diffuser with a proposed area ratio of 16 to 1 appeared 
highly problematical from a recovery point of view. In addition, the effects on 
the flow of an instrument -mounting employed inside the parallel pipe appears 
theoretically unpredictable. Finally, all of these effects combined presents the 
problem of matching the pressure changes and flow rates required from the system 
against the actual quantities being available by the induction under operational 
conditions. 

TEST EQUIIMENT 

The test equipment essentially consiste^^of a parallel pipe which was provided 

% 

with a well-rounded entry at its inlet and was ''onnected to a diffuser at its outlet. 


* The actual scaling factor was 1:9.6. 
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Located upstream from the entry a venturi meter of standard design was employed 
for measuring the flow rate. The general arrangement is shown in figure la. 

The parallel pipe was made of extruded acrylic tubing with approximately 
.00476 m wall thickness and was fitted with flanges at each end. To establish 
the exact location of the section where diffusion started, the end of the pipe 
was slightly tapered from about one inch distance from its exit plane. Pressure 
tappings were distributed at regular spacing along the pipe but the distance 
between tappings was decreased towards the ^ipe exit where large pressure gradients 
were anticipated. Some sections were provided with a single tapping while others 
had three taps which were interconnected, thus allowing for the average of three 
pressures to be recorded in any one plane of measurement. 

While there was only one pipe employed throughout the experiments, there were 
two diffusers tested. The first diffuser was made of plastic material and had 
straight walls with a taper angle of 3 degrees, while the second diffuser was made 
of sheet metal. This started with a lower taper angle of 2.4 degrees which changed 
to a higher angle of 4.76 degrees halfway along the diffuser length, as shown in 
figure lb. Pressure tappings were distributed along the diffusers at regular 
spacing. 

Tests were conducted under two different operational modes: (a) under suction, 

and (b) under pressure, and the test arrangement was changed accordingly. Under 
suction , the flow was "inducted" and operation was maintained with the aid of a 
vacuum system consisting of an air ejector to which the primary air was supplied 
jy a compressor. The ejector was specially designed for those tests to reduce pres- 
sure to about 22,86 cm mercury below atmosphere inside two tanks each of 28.317 m^ 
capacity. Evacuation time to the desired vacuum lasted for about 15 minutes. 

Under suction operation all pressures along the tube and diffuser were measured 
with a multitube manometer, while the suction pressure in the vacuum line was read 
on a single tube manometer, both manometers using mercury as indicating fluid 
(see fig. 2a). Under pressure , flow through the system was maintained by a 
powerful centrifugal blower which discharged into a settling chamber where the 
pressure could be varied from .0508 m to .254 m water gage. Fluctuations of the 
airstream were damped out with honeycomb screens. From the side of the settling 
chamber the flow entered the parallel pipe through the inlet section and was 
finally discharged into the atmosphere through the diffuser (see fig. 2b). 

Under the pressure operation, velocity traverses across various cross sections of 
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the pipe and the diffusers were obtained with a standard pitot static tube of 
,001587 m diameter. 

For velocity traverses and for measuring the flow rate, inclined tube mano- 
meters were employed using alcohol as indicating fluid. Atmospheric temperature 
and pressure were read from a standard barometer located near the test equipment. 

For the boiindary layer suction experiments in the second diffuser a small 
blower unit was used and attempts at sucking the layer were made both at diffuser 
inlet and at the halfway section. The flow rate was measured with a small Venturi 
meter built into the suction line. 


PROCEDURE 


The experiments were performed under either quasi-steady or fully steady 
flow conditions. More particularly, under suction the induced aiiflow into the 
tanks through the pipe was greater than the evacuating volume flow, because the 
compressor was unable to maintain an adequate supply of primary air to the ejector 
under decreasing vacuum. However, when the flow attained sonic speed at the pipe 
exit the pressure distribution along the pipe remained unchanged for some short 
period of time, long enough for it to be considered steady, while under subsonic 
flow conditions at pipe exit the pressures were steadily rising. In either case, 
the pressure distribution was instantaneously recorded by photographing the mano- 
meter board. On the other hand, fully steady operation was maintained when 
operating under the pressure mode. 

For the suction mode experiments the following routine applied. Air was 
withdrawn from the closed tank until a maximum vacuum of about 23,368 cm Hg was 
attained. Each test was started by quickly opening the admission valve to the 
tank 'ind immediately observing the manometer board. When the desired flow condi- 
tions were attained the board was photographed and simultaneously the suction 
pressure and the pressure differential of the Venturi were recorded. Tlie tests 
were subsequently repeated for several different suction pressures. Immediately 
after each test the admission valve was closed and remained closed until maximum 
vacuum in the tank was restored. 

Several test runs were performed with unrestricted flow through the pipe and 
one run was performed with a protruberance that restricted the passage at a parti- 
cular cross section. It was located .2794 m from the |)ipc entrance. The protruber- 
ance was shaped after an instrument holding device and it simulated partial blockage. 
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During the pressure mode experiments, velocity traverses were made at the 
relevant points along the flow — at the inlet to the parallel pipe, at the inlet 
to the diffuser, at the half-way section of the diffuser, and at the exit. These 
tests were subsequently repeated with boundary layer suction which was applied 
to either the inlet to the diffuser or the half-way section of the diffuser. 
During these experiments the effects resulting from varying the suction flow 
rate were also studied. During the tests the atmospheric temperature and 
pressure were recorded. 


EVALUATION OF TEST RESULTS 

For calculating the Mach number M from the mass flow rate, stagnation 
temperature and static pressure, the Fanno equation was used 



Details of the procedure for calculating M arc given in reference 1. 

For the calculation of the recovery in the diffuser, the recommended 
expression was employed (ref. 2). 


c 

P 


P2 " Pi 

J PlVl^ 


( 2 ) 


where Vj = /KgRDj and = T^/ ^1 + j value of Pj can be 

found from the equation of state pj = Pi/gRB^. 

It is noted that the entry condition to the diffuser (subscript 1) was 
assumed to be identical to the exit condition from the parallel pipe. 

For calculating the loca'* skin friction along the parallel pipe it was assumed 
that the distances between consecutive pressure tappings were small enough for the 
factor f to be considered constant over those distances. Then 


f = 


D 

4AL 




K + 1 
2K 


lu 




( 3 ) 


b 



where 


1 + 


Ml 


M,2 


and 


4 ) = 

M2 M22 

Here the subscripts 1 and 2 refer to the inlet and outlet section respectively 
of the pipe length (AL) under consideration. 

The mass flow rate of air through the Venturi was calculated with the 
expression 

W = Cj a^ /ZgyAp (4) 

where the pressure differential AP was measu»'cd on the inclined manometer 
connected to the throat section. The other limb of the manometer was open to 
the atmosphere. 


DISCUSSION OF THF: RESULTS 

The results of the experiments are presented in figures 3 to 8. In figure 3, 
the distribution of Mach number along the pipe and diffuser is shown for four 
different suction pressures from 12.7 cm to 21.1 cm Hg vacuum. In figure 4 the 
pressure recovery of the test diffusers is plotted against inlet Mach number, 
while in figure 5 the variation of the friction factor along the pipe is shown 
for the specific case when M = 1 at exit. In figures b and 7 the effects of 
boundary layer suction at various locations arc shown and in figure 8 the velocity 
distribution across the inlet plane to the parallel pipe is presented. 

Figure 3a shows that the Mach numlier at inlet to the parallel pipe was about 
0.62 as against 0.58 calculated with an assumed constant friction factor of 0.02. 
About halfway downstream from inlet the flow attained M = 0.7, and from x/d = 20 
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onwards there appears a rapid increase in M, The gradient dM/dx is highest at 
the exit, where M becomes just slightly greater than unity*. In going downstream 
a rapid decrease in M is experienced in the inlet region of the diffuser and ulti- 
mately the flow exits with M “ 0.05. For suction pressures lower than 18 cm Hg 
the pressure distributions change, resulting in lower Mach numbers as shown in fig- 
ures 3c and 3d. While at pipe exit the Mach number falls substantially below unity 
the curves remain rather similar in appearance. After the insertion of the piotrub- 
erance some tests were repeated to study the effect of partial blockage. The 
changes were found negligible and presentation of these results has been omitted. 

The pressure recovery factor for the original diffuser was found to be 

about 0.8 (fig. 4) and was slightly improving with increasing inlet Mach number. 

This figure is in good agreement with results obtained by other investigators 
(refs. 3 and 4). The pressure recovery for the proposed new diffuser was found 
about 4 percent higher and the experiments show that most of the recovery is 
already accomplished at halfway along the diffuser. 

The local friction factor varies considerably along the pipe (fig. 5). It 
decreases first, then increases and finally markedly decreases near the exit. The 
scatter of the experimental points suggests that under unsteady flow the pressure 
measurements obtained by the mercury-in-glass piezometers have rather limited 
accuracy. It becomes clear, however, that the assumption of constant friction 
factor in a highly accelerating gas flow docs not hold. 

Velocity traverses obtained at various cross sc tions, namely at the inlot 
(X), outlet (ZJ, and mid-section (Y) of the second half of the new diffuser show 
the beneficial effects of boundary layer control with suction. The suction was 
applied at mid-section X only and sliows little I if any) effect on the profile 
at X (fig. ba) . The effect, however, is more noticeable at section Y (fig. bb) 
and is quite marked at the exit section Z. The gradual improvement of profiles 
with increasing boundary layer suction is shown in figure 7. 

Finally, the velocity distribution obtained at the inlet to the parallel pipe 
(fig. 8) shows the effectiveness of the contraction. While the velocity distribution 
is fairly uniform across the center portion of the pipe, it appears to be slightly 
higher on both sides due to curvature of the contraction walls. This effect is 


The increase of M alcove unity at the exit plane is attributed to the streamlines 
slightly diverging over a short distance aliead of the inlet to the diffuser. 
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considered negligible and usually disappears shortly after the flow enters the 
parallel pipe. 

From the measured local friction factor an overall can be established and with 
this value assumed constant along the pipe, a variation of the inlet with the exit 
Mach number may be predicted as shown in figure 9. It appears that at lower Mach 
numbers (say, up to 0.4) the flow may be considered as incompressible. 

coNci.nsroN 

The general conclusion tliat may be drawn from the experiments appears favor- 
able and the projected duct design should meet the requirements. In other words, 
by scaling up the experimental results it may be predicted that the operational 
demand from the duct appears to match the rerformance of the induction system. 
Relevant details of calculations arc prcsep':»‘d in Appendix A. 

More particularly it was found that: 

1. A Mach number equal to unity at the exit from the parallel pipe may be 
anticipated in the range of induction pressures 18.03 to 21.1 cm Hg (vacuum) corres- 
ponding to 7888 kgf/m^ and 7473.6 kgf/m^, respectively. Since the effect of stream 
line curvature on the exit Mach number is rather difficult to estimate accurately, 

a conservative value that may be recommended for the induction pressure would be 
about 20.32 cm Hg corresponding to about 7593 kgf/m^. 

2. The pressure recovery in both of the diffusers may be considered satis- 
factory. However, the velocity distribution at the exit of the new proposed 
diffuser appears more attractive especially when boundary layer suction is 
applied at half distance along its length. Row oscillations, whicli always 
accompany separated flow, thus ma> be reduced or evi-n eliminated because the 
risk of flow separation is minimized with "full" velocity profiles. Relevant 
details of boundary layer suction are presented in Appendix 8. 

3. Variation of the friction factor along the pipe is ctnisiderable and as 
a result, the experimentally found Mach number at pipe inlet was found liigher 
than the calculated figure based on constant friction. 

4. The design of the contraction (air intake) appears satisfactory as the 
velocity distribution across the pipe inlet is fairly uniform. 




5, A comparison drawn between pressure distributions along the pipe with 
unrestricted flow and restricted passage due to the protruberance indicated no 
substantial difference. Thus the partial blockage, due to the instrument panel, 
may be predicted to be negligible. 
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(a) Parallel pipe (all dimensions in meters), 
l igure 1. (icneral Arrangement of Test Duct-systems . 




SECOND, PROPOSED NEW DIFFUSER 


showing location of pressure tappings (all dimensions in meters) 
Figure 1. Concluded. 







(a) Set-up for suction mode operation. 
Figure 2. Photograph of experimental apparatus. 
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BOUNDARY LAYER SUCTION 
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(b) Set-up for pressure mode operation 


Figure 2. Concluded 


pFPT?OTlTTrrf.;iT,firf 0 







(a) Suction pressure .211 m Hg. 


Figure 3. 


Variation of Mach Number along Parallel Pipe. 





(b) Suction pressure .180 m Ilg. 


Figure 3. Continued. 




Length Along Pipe and Diffuser, x/d 

(c) Suction pressure .152 m Hg. 


Figure 3. Continued, 




Mach Number 



(d) Suction pressure .127 m Hg. 


Figure 3, Concluded. 



First half of 
Proposed Diffuser 



Inlet Mach Number, M 


(b) Proposed new diffuser. 


Figure 4. Concluded. 




Length Along Parallel 


Figure 5, Variation of friction factor along pa 
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Figure 6. Variation of velocity across sections of diffuser. 
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(b) Second diffuser. Section Y. 
Figure 6. Continued. 
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(c) Exit Plane, Section Z. 
Figure 6. Concluded. 
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Figure 7. Variation of velocity across exit plane (Z) for different suctions. 




Figure 8. Variation of velocity across inlet plane to the 
parallel pipe. 
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APPENDIX A 


CALCULATION OF THE DESIGN POINT FOR COMPRESSOR 


From the experiments performed on the model, the design operational condition 
representing the anticipated mass flow rate and diffuser exit-pressures can be 
estimated for the full scale duct. 

Assuming isentropic flow, the stagnation pressure at inlet to the parallel 
pipe would be atmospheric, with the experimentally found = 0,62 at inlet, 
the static pressure 

^ ^o (14, 73(6894.75) 

(l + Mi2 j (9.80665) [l + Mj2 j 


(14.7) (6894.75) 


(9.80665) 1 + .2(.62)2 




= 7975.25 kgf/m2 


Since the flow is frictional and M 2 = 1 


at pipe exit, from the Fanno equation 



hence 


P2M2 



P2 = 



= 4686 kgf/m2 


The static temperature at pipe exit at t^ = 70® F, 




By the equation of state p/y = R9, 


P2 


- m 
kg ”R 


(53.3) (4.448) (.3048J 
(9.80665) (.45559) 


(4686) (9.80665) (.45359) 

^2 ^ (53.3)(4.448\.3048)(441.7) 


Y 2 = 0.653 kg/m^ . 


Velocity of air at sonic condition for y = 1,4, 

_ (55.5) (4.448) (.3048) _ kgf - m 

(9.80665) (.45359) “ kg ”R 

g = (32.2) (.3048) = 9.81456 m/scc^ 

V 2 = ^(1.4) (9.81456) (16.245)0 = 314.02 m/sec 

Hence the mass flow rate, with pipe area 
A = |^(1.25)(.0254)j^ J = .007917 
and without taking the boundary layer into consideration 
W = yAV = (.653) (.007917) (314.02) = 0.16254 kg/scc. 
Experimentally, mass flow rate of the model found 
(.3515) (.45359) = .159437 kg/sec . 


Without the boundary layer thickness the sealed up flow rate 
would be W = (.653) (.785) (.6j^.903) (314.2) = 14.95 kg/sec; with a 
layer displacement, this reduces to an actual flow rate 


in the duct 
scaled boundary 


W =14,817 kg/scc . 
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Since sonic flow was attained at about ,203 in mercury, the exit pressure from the 
diffuser was found 


P 


(14,7) (6894,75) 
9,80665 


= 7579,08 kgf/m2 . 


The point of operation is shown on the compressor performance chart (fig, Al) 
as OP, 
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, discharge pressure, kgf/m^ 




APPENDIX B 


For the estimation of boundary layer suction power the following approximate 
calculation was employed. 

For the maximum volume suction for the model AQ = (. IDS) (,02851hS) = 

,00297 m^/sec there corresponds a AQ = (92. 16) (.00297) = .274016 mV-sec in 
the duct. Experiments show that the pressure at halfway section was 7733.76 kgf/ 
hence the difference between atmospheric pressure Ap = 10335.11 - 7733,76 = 
2601.35 kgf/m2. 

The theoretical power required to remove this air is 


p = (A4>)(AP) = (.274016) (2601.35) = 712.81 kgf - m/sec = 6.995 kW. 


To match this requirement a 11.60 kW motor would be recommended. 


